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Using the recent WMAP determination of the baryon-to-photon ratio, 10 10 i] = 6.14 to within 
a few percent, big bang nucleosynthesis (BBN) calculations can make relatively accurate 
predictions of the abundances of the light element isotopes which can be tested against 
observational abundance determinations. At this value of 77, the 7 Li abundance is predicted 
to be significantly higher than that observed in low metallicity halo dwarf stars. Among the 
possible resolutions to this discrepancy are 1) 7 Li depletion in the atmosphere of stars; 2) 
systematic errors originating from the choice of stellar parameters - most notably the surface 
temperature; and 3) systematic errors in the nuclear cross sections used in the nucleosynthesis 
calculations. Here, we explore the last possibility, and focus on possible systematic errors 
in the 3 He(a,7) 7 Be reaction, which is the only important 7 Li production channel in BBN. 
The absolute value of the cross section for this key reaction is known relatively poorly both 
experimentally and theoretically. The agreement between the standard solar model and 
solar neutrino data thus provides additional constraints on variations in the cross section 
(5*34). Using the standard solar model of Bahcall, and recent solar neutrino data, we can 
exclude systematic S34 variations of the magnitude needed to resolve the BBN 7 Li problem 
at > 95% CL. Additional laboratory data on 3 He(a,7) 7 Be will sharpen our understanding 
of both BBN and solar neutrinos, particularly if care is taken in determining the absolute 
cross section and its uncertainties. Nevertheless, it already seems that this "nuclear fix" to 
the 7 Li BBN problem is unlikely; other possible solutions are briefly discussed. 



1 Introduction 



The recent all-sky, high-precision measurement of microwave background anisotropies by 
WMAP has opened the possibility for new precision analyses of big bang nucleosynthesis 
(BBN). Until now, one could use the predictions of standard BBN j2JIH] for the abundances of 
the light element isotopes, D, 3 He, 4 He, and 7 Li and compare those results with the observa- 
tional determination of those isotopes and test the concordance of the theory. If concordance 
is found, the theory is also able to predict the value of the baryon-to-photon ratio, rj. Indeed, 
concordance is found, so long as a liberal estimation of systematic uncertainties are included 
in the analysis. The accuracy of the predicted value of rj from BBN alone based on likelihood 
methods (3 El El Ej is modest: i]i = 5.7tl'^ when D, 4 He, and 7 Li are used, and 7710 = 6-0-o.i 
when using D alone, where rjxo = 10 10 r]. This pales in comparison with the recent WMAP 
result of Vt B h 2 = 0.0224 ± 0.0009 which is equivalent to r/i , C MB = 6.14 ± 0.25. This result 
is the WMAP best fit assuming a varying spectral index and is sensitive mostly to WMAP 
alone (primarily the first and second acoustic peaks) but does include CBI [8 and ACBAR 
jHj data on smaller angular scales, and Lyman a forest data (and 2dF redshift survey data 
[TU] ) on large angular scales. 

If we use the WMAP data to fix the baryon density, we can make quite accurate predic- 
tions for the light element abundances. At this density, we can make a direct comparison 
[TT] between theory and observation as shown in Table ^ 

Table 1: Light Element Abundances: BBN Predictions and Observations 



element 


theory 


Observation 


D/H 
4 He 
7 Li 


2.75i^ 4 x lO" 5 

n 9_8_+ - 0004 
U.2484_ . 00 o 5 

3.82ljj$ x lO" 10 


2.78 ±0.29 x 10" 5 
0.238 ±0.002 ±0.005 
1.23 + °'? 4 x lO" 10 



As one can see, the agreement between the predicted abundance of D/H and the observed 
value (based on the average of the 5 best determined quasar absorption system abundances 
[121 Q__ EEH EH] ) is perfect. The comparison with 4 He is less good, as BBN predicts a mass 
fraction which is high compared to most observations fOJ El __H1 EI] ■ The value in Table Q 
is based on a combined analysis ^H] which is close agreement with the recent observations of 
|19j . One should note that 1) the data of alone give a higher value for the 4 He abundance 
Y p = 0.242 ± 0.002 ± 0.005, and 2) important systematic effects have been underestimated 
20 . Among the most probable cause for a serious underestimate of the 4 He abundance 
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is underlying stellar absorption. Whether or not this effect can account for the serious 
discrepancy now uncovered remains to be seen. 

Clearly the key problem concerning the concordance of BBN theory and the observational 
determinations of the light element abundances is 7 Li. The primordial abundance of 7 Li is 
determined from the "Spite plateau" [21] in Li/H observed in low metallicity halo dwarf stars 
(extreme Population II). The observed value is clearly discrepant with the BBN+WMAP 
prediction. The cause of the discrepancy may be: 

• Stellar depletion of 7 Li - however, the lack of dispersion in the observed data, make it 
unlikely that dispersion alone can account for the difference. 

• Stellar parameters - the determined 7 Li is sensitive to the assumed surface temperature 
of the star. However, to account for a discrepancy this large temperatures would have 
to be off by at least 500K. This may not be reasonable. 

• The nuclear rates - this is the case we wish to explore here. 

Of course, it is also possible that the 7 Li discrepancy is real, and points to new physics. 
However, it is our view that at present, the case for new physics is not compelling, though 
it certainly merits serious investigation. Furthermore, a firm rejection of the more "prosaic" 
possibilities we have outlined is a prerequisite which must be satisfied before we are driven 
to more radical and exciting new solutions. It is in this spirit that we investigate possible 
systematic errors in the BBN theoretical predictions for 7 Li. 

Uncertainties in the nuclear reaction rates which determine 7 Li are dominated by 3 He(a, 7) 7 Be. 
There has been concerted experimental and theoretical effort to understand this reaction, 
and indeed the cross section shape versus energy appears to be well-understood [22]- How- 
ever, a challenge to both experimental and theoretical work has been the determination 
of the absolute normalization of of the cross section. This uncertainty propagates into an 
overall systematic error in the 3 He(a,7) 7 Be rate. 

We thus pose the following question. Independent of the quoted (or derived) laboratory 
uncertainties in 3 He(a,7) 7 Be, what is the maximum allowable amount that this rate can be 
adjusted. Of course, we are not completely free to adjust this rate, since this nuclear reaction 
occurs in the Sun and is in part responsible for the observed flux of solar neutrinos. Thus our 
goal is to use the standard solar model [22] as a constraint on the BBN nucleosynthesis rates. 
In order to reduce the predicted 7 Li abundance in Table 1, to the observed one requires a 
reduction in the 3 He(a,7) 7 Be by a factor of 0.27. We show that by using the concordance 
between the standard solar model and the observed flux of solar neutrino, this is excluded 
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at the 99.9999 % CL. At the 95% CL, the largest reduction factor possible is 0.65. Thus, it 
is not possible to argue that the uncertainties in nuclear reactions are solely responsible for 
the 7 Li discrepancy. 

In section 2, we detail the problem of BBN produced 7 Li. In section 3, we discuss the key 
nuclear reactions which contribute to the overproduction of 7 Li. We derive our constraints 
on these reactions using the observed flux of solar neutrinos in section 4. A summary and 
discussion is given in section 5. 

2 The Overproduction of 7 Li 

As noted in Table 1, the BBN 7 Li abundance is predicted to be 3.82^o;<3q x 10~ 10 for 7710 = 
6.14 ± 0.25. This result [TT] is based on a BBN calculation [7J using the updated rates 
compiled by the NACRE collaboration [23]. Other calculations tend to give even higher 
values, e.g., 7 Li/H = 4.87ig;{£ j2EJ; 7 Li/H = 4.18 ± 0.46 [27]. These results differ due to the 
different nuclear data sets and procedures used to fit them and derive thermonuclear rates. 
The variations are thus a measure of known systematics in the 7 Li prediction. 

The observed Li/H value in Table 1 reflects the inferred mean abundance in the atmo- 
spheres for a set of Pop II stars. The analysis is that of [2H|, based on the data of j2H]- The 
data sample consists of 23 very metal poor halo stars, with metallicities ranging from [Fe/H] 
= -2.1 to -3.3. The data show a remarkably uniform abundance of Li and negligible disper- 
sion about a tiny slope which is consistent with the production of some Li in Galactic cosmic 
ray collisions (primarily a + a). Note that any Galactic component of Li only compounds 
the BBN discrepancy. 

The 7 Li value in Table 1 assumes that the Li abundance in the stellar sample reflects 
the initial abundance at the birth of the star; however, an important source of systematic 
uncertainty comes from the possible depletion of Li over the > 10 Gyr [HO] age of the Pop 
II stars. Stellar interiors can burn Li and alter its surface abundance. The atmospheric Li 
abundance will suffer depletion if the outer layers of the stars have been transported deep 
enough into the interior, and/or mixed with material from the hot interior; this may occur 
due to convection, rotational mixing, or diffusion. However, if mixing processes are not 
efficient, then Li can remain intact and undepleted in a thin outer layer of the atmosphere, 
which contains a few percent of the star's mass but is the portion of the star's material that 
is observable. 

Standard stellar evolution models predict Li depletion factors which are very small 
(<0.05 dex) in very metal-poor turnoff stars [3T]. However, there is no reason to believe 
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that such simple models incorporate all effects which lead to depletion such as rotationally- 
induced mixing and/or diffusion. Current estimates for possible depletion factors are in 
the range ~ 0.2-0.4 dex jH2]- While the upper end of this range is close to the required 
depletion factor of ~ 0.3 necessary to account for the difference in the BBN and observed 
abundance, depletion models typically predict the existence of star-to-star differences in ob- 
served Li abundances due to the range of stellar rotation and other intrinsic stellar properties 
to which the models have some sensitivity. As noted above, this data sample [21] shows a 
negligible intrinsic spread in Li leading to the conclusion that depletion in these stars is as 
low as 0.1 dex. 

Another important source for potential systematic uncertainty stems from the fact that 
the Li abundance is not directly observed but rather, inferred from an absorption line strength 
and a model stellar atmosphere. Its determination depends on a set of physical parameters 
and a model-dependent analysis of a stellar spectrum. Among these parameters, are the 
metallicity characterized by the iron abundance (though this is a small effect), the surface 
gravity which for hot stars can lead to an underestimate of up to 0.09 dex if log g is overes- 
timated by 0.5, though this effect is negligible in cooler stars. Typical uncertainties in log g 
are ±0.1 — 0.3. The most important source for error is the surface temperature. Effective- 
temperature calibrations for stellar atmospheres can differ by up to 150-200 K, with higher 
temperatures resulting in estimated Li abundances which are higher by ~ 0.08 dex per 100 K. 
Thus accounting for a difference of 0.5 dex between BBN and the observations, would require 
a serious offset of the stellar parameters. 

We note however, that a recent study [33] with temperatures based on Ha lines (consid- 
ered to give systematically high temperatures) yields 7 Li/H = (2.19 ± 0.28) x 10~ 10 . These 
results are based on a globular cluster sample and do show considerable dispersion. A related 
study (also of globular cluster stars) gives 7 Li/H = 2.29 x 10 -10 [HI]. The difference between 
these results and the BBN value is just over 0.2 dex making it plausible that depletion may 
be responsible for the difference in these stars which show systematically high temperatures. 
It remains an open question why stars in a globular cluster-which are usually thought of 
as sharing a common origin site and epoch-seem to show a larger Li dispersion (and higher 
temperatures) than field halo stars whose evolution has not been so tightly related. 

Finally, the remaining source of systematic uncertainty pertains not to the observations, 
but to the BBN calculation itself. Here we will limit ourselves to a discussion of those cross 
sections which have a bearing on the production of 7 Be, which is the dominant source of 
mass-7 at the high values of i] consistent with the WMAP result. 1 As such the principle 

1 At ijio = 6.14, the production ratio is 7 Be/ 7 Li = 10.8. Of course, the 7 Be eventually suffer electron 
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cross sections of interest are: 3 He(a,7) 7 Be and 7 Be(p, 7) 8 B. The reaction 7 Be(ra,p) 7 Li is not 
of interest since it does not largely affect the final abundance of 7 Li. 

3 Nuclear Rates contributing to BBN 7 Li production 
3.1 Standard BBN 

Since our aim is to fix the 7 Li problem by changing nuclear reaction rates, specifically the 
3 He(a, 7) 7 Be and 7 Be(p, 7) 8 B reactions, it is important to understand how they do or do not 
impact primordial nucleosynthesis. We will start with the all-too-familiar n(p, 7)0? reaction 
and how it affects the light element yields. This will guide us when looking specifically at the 
other reactions. It is well-known that nucleosynthesis in the early Universe is delayed due 
to the deuterium bottleneck. It is important to understand how the deuterium bottleneck 
affects the abundances of the light elements. The delay being caused by the large number of 
photons to baryons, which makes the deuterium photo-destruction rates much larger than 
the production rates. At lower temperatures, about 70 keV, deuterium production proceeds 
and the burning into heavier nuclei occurs until the Coulomb barrier halts nucleosynthesis. 
We burn until we deplete the neutron fuel and the Coulomb barrier stops charged-induced 
reactions, happening at a temperature around 50 keV. 

While the bottleneck is in place, neutrons and protons remain at their weak freeze-out 
values, except for the occasional n-decay, and deuterium at its equilibrium value. The other 
light element abundances exist in a quasi-static statistical equilibrium, being determined by 
various algebraic combinations of the important thermonuclear reaction rates [35J. When 
the bottleneck ends and the neutron fuel is depleted into 4 He, these abundances tend to 
freeze-out at particular values depending on the overall baryon content in the universe and 
when the bottleneck ended. 

The deuterium bottleneck ends when the photo-dissociation rate ^(7, n)p becomes less 
significant than the rap-capture rate. This is done by setting the ratio of equilibrium abun- 
dances to unity, X d /X p X n ~ 1 for which there is an analytic expression [35] that can be 
solved iteratively: 

rp _ ^± f]\ 

d 28.7 - In (7710/6.O) - 1.5 In (T d /MeV) 1 ; 

One finds that for = 2.224 MeV and 1] W = 6.0, the deuterium bottleneck ends at = 0.07 
MeV. This method can also be used to determine when certain processes become inefficient, 
capture and decay to 7 Li before recombination and long before incorporation into Pop II stars. 
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such as 7 Be(7,a) 3 He and 8 B(7,p) 7 Be. Using equations that can be similarly solved itera- 
tively, we find the temperature for which the photo-dissociation of 7 Be and 8 B listed earlier 
is not significant: 

34 28.2 - In (77x0/6.0) - 1.5 In (T 34 /MeV) 1 ; 

17 30.1 - In (7710/6.O) - 1.5 In (T 17 /MeV) 1 ; 

Using Q- values of 1.587 MeV and 0.137 MeV for the 7 Be(7,a) 3 He and 8 B(7,p) 7 Be respec- 
tively, we find T34 = 0.05 MeV and T 17 = 0.004 MeV. These temperatures are below the 
deuterium bottleneck, suggesting that if these reactions dominate the formation of these 
elements, then they should exist at their equilibrium abundances. This is true for the 
3 He(a,7) 7 Be reaction. However, 7 Be(j9, 7) 8 B must compete with 8 B decay, with mean life- 
time ~ 1 s. An equilibrium abundance of 8 B during the epoch of big bang nucleosynthesis 
will be greatly suppressed and in fact is completely negligible. In other words, 8 B and its 
decay products do not contribute to the primordial light element abundances. Knowing that 
3 He(a,7) 7 Be is the dominant contribution to the BBN 7 Li abundance prediction, we now 
discuss what we need to fix the 7 Li problem. 

3.2 Modified BBN: A Nuclear Solution to the Lithium Problem 

The question of interest to us here, is to what extent can these key rates be altered to enhance 
the 7 Be ( 7 Li) abundance and yet remain consistent with experimental constraints. To this 
end, we define a new S-factor 2 Si 7 EW which we assume for simplicity to be proportional 
to the old one, S® 7 LD . Note that for Sn, a proportionality factor between and 2 does 
not change the BBN predictions significantly. In contrast, the dependence of the mass 7 
abundance on S34 is nearly 1:1, as apparent in Table El and in good agreement with the 
results of [6 , who find that 7 LiBBN oc S® 4 95 . 

As discussed above, there are two sets of 7 Li observations we can try to match by renor- 
malizing the 3 He(a,7) 7 Be reaction. Using the 7 Li measurements of a metal poor globular 
cluster [HI] would require a change of S^ 4 EW = 0.53S^ 4 LD . Using the 7 Li measurements of 
metal poor stars in the Galactic halo j2Hj would require a change of S^ A EW = 0.27S^ 4 LD . 

The determination of the BBN light element yields is from |7 , where new normalizations 

and errors to the NACRE [23] rates important for primordial nucleosynthesis have been 

assigned. The value, S? 7 LD (0) = 0.021 ± 0.002 keV b is taken straight from the NACRE 

2 The S'-factor is defined by the cross section: S(E) = a(E)E exp(87r 2 aZ 1 Z 2 /w). The last term is the 
Coulomb penetration factor, in which Zi are the charges of the incoming nuclei and v their relative velocity. 
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Table 2: 7 Li Sensitivity to S34 



aNEW 1 oOLD 

^34 / °34 


yNEW jyOLD 


1.00 


1.00 


0.75 


0.76 


0.50 


0.52 


0.25 


0.28 



collaboration. For 3 He(a,7) 7 Be, the BBN calculation uses the renormalized NACRE rate 
S^ LD (0) = 0.504±0.0534 keV b 7\. As one can see, shifts in the 3 He(a, 7) 7 Be cross section as 
large as that necessary to produce S^ w are strongly excluded given the cited uncertainties 
for this reaction. Although adjustments in the nuclear cross-sections of this size are unlikely 
given the stated experimental errors, one could worry that additional systematic effects are 
present, particularly given the difficulties in establishing the absolute normalization for this 
reaction. As stated in the Introduction, these rates in particular can be bounded by another 
means. In the next section, we will determine the maximum possible downward adjustment 
to S34 which is consistent with solar neutrino fluxes. 

The effect of changing the yields of certain BBN reactions was recently considered by 
Coc et al. [2Z1- I n particular, they concentrated on the set of cross sections which affect 7 Li 
and are poorly determined both experimentally and theoretically. In many cases however, 
the required change in cross section far exceeded any reasonable uncertainty. Nevertheless, it 
may be possible that certain cross sections have been poorly determined. In j27| , it was found 
for example, that an increase of the 7 Li(cZ, n)2 4 He reaction by a factor of 100 would reduce 
the 7 Li abundance by a factor of about 3. Another reaction which is poorly determined 
is 7 Be(<i,p)2 4 He. An increase in this rate by a factor of ~ 100 could also alleviate the 7 Li 
discrepancy 

4 The Sun as a Nuclear Laboratory 

The 3 He(a, 7) 7 Be reaction plays a crucial role not only in BBN 7 Li synthesis, but also in solar 
neutrino production. In particular, this reaction is responsible for the creation of 7 Be, which 
will then either (1) produce a monoenergetic neutrino via electron capture 7 Be(e~, z/ e ) 7 Li, or 
(2) produce 8 B via radiative capture of a proton, 7 Be(p, 7) 8 B. The branching between these 
paths determines the solar 8 B abundance and thus directly sets the flux of 8 B neutrinos. 
SNO (as well as Super-K) are sensitive exclusively to these neutrinos. Furthermore, SNO 
measures directly the total 8 B neutrino flux, with no assumptions about mixing j^B]. They 
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find: 

8 = [5.21 ± 0.27(stat) ± 0.38(syst)] x 10 6 cm" 2 s _1 , (4) 

where this is determined with no assumed shape of the 8 B energy spectrum. This flux thus 
offers a constraint on the 3 He(a,7) 7 Be reaction, as follows. 

The Standard Solar Model of Bahcall |3Ij can be used to predict the solar neutrino fluxes 
that can be observed by experiments. These fluxes depend upon various solar parameters, 
such as the luminosity, the chemical abundances, and nuclear fusion cross sections. In fact, 
the neutrino flux uncertainties are dominated by the cross section errors. Provided by j37| . 
simple scalings between neutrino fluxes and these cross sections robustly describe the SSM 
predictions. The 8 B neutrino flux scaling is: 

8 oc S^S^^SgSj*. (5) 

Here, the S"s are the astrophysical S-factors, except for S e7 . The S e7 reaction is the electron 
capture rate on 7 Be. One usually takes some nuclear rate compilation and uses the ^-factors 
to evaluate the neutrino flux given these scalings. Two such nuclear compilations are from 
Adelberger et al [21] an d the NACRE collaboration [23]. Their determinations relevant for 
this work are shown in Table El 



Table 3: Shown are the results from the nuclear fusion rate compilations. 



Reaction 


Adelberger [24J (keV b) 


NACRE [25 j (keV b) 


p(p, e + z/ e ) 2 H 


S n = 4.0 x 10- 22 (1.0 ± 0.007±^ 2 ?) 


S u = 3.94 x 10" 22 (1.0 ± 0.05) 


3 He( 3 He, 2p) 4 He 


S 33 = 5.4 x 10 3 (1.0 ± 0.074) 


S 33 = 5.18 x 10 3 (1.0 ±0.06) 


3 He(a,7) 7 Be 


S u = 0.53(1.0 ±0.09434) 


S 3A = 0.54(1.0 ±0.167) 


7 Be(jo,7) 8 B 


S u = 0.019(1.0l£?i) 


S 17 = 0.021(1.0 ±0.11) 


7 Be(e-,^ e ) 7 Li 


S e7 = 5.6 x 10- 9 (1.0±0.02) s" 1 


-NA- 



Our strategy is thus to use the SNO measurements of the 8 B neutrino flux and the SSM 
to constrain S34. This is accomplished via the scalings in eq. (JHJ). A complication is that 
these scalings also depend on other reactions, none of which are significant for BBN, and 
all of which are better measured than 3 He(a,7) 7 Be. This approach amounts to the extreme 
case in which we ignore all of the hard- won laboratory and theoretical information on £34, 
using only solar neutrino data as well as constraints on other reactions, Sn. This can be 
viewed as providing independent information about S34, or as a test of the systematics in 
the normalization, which is a salient feature for the BBN 7 Li problem. Our results will thus 
use the Sun to provide new and independent limits on the systematics of Su- We will derive 
these using both approximate analytical methods and more accurate numerical methods. 
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4.1 Analytic Formalism and Results 

We can estimate the impact these rate compilations have on the neutrino flux, by doing 
linear error propagation as follows: 

We find that the Adelberger and NACRE compilations predict a 8 /0 8 = ±0.19 and cr 8 /08 = 
±0.22 respectively using this linear approximation. With these results, we find that the error 
in the predicted flux is determined primarily by the Sn, S34 and Su reactions. With our 
ultimate aim of constraining S34, we will have to treat at least the Sn and Su uncertainties 
directly, in addition to the error in the solar neutrino flux measurement. 

We can now use the scalings in eqn.|5]to estimate the likely value of S34 and its uncertainty, 
based on the SNO observations 

s* = (Su\ i21 fs^\ 0A9 fs^Y 1 - 23 fs^y* 

£34,0 \S\ifi ) \S33fl J \Sn 7 o J \S e 7fiJ \4>8,o 

where we use the Bahcall et al. results for the Adelberger and NACRE reaction complications 
(the Si t o) to determine the flux normalization 

</,adl = 5m x 10 -6 cm -2 s -i ( 8 ) 

0NAC = 5 44 x 10 -6 cm -2 s -l (9) 

These normalizations are both in excellent agreement with the observed flux (eq. EJ). In the 
extreme case in which all of the small mismatch between predicted and observed fluxes is 
attributed to S34, we expect a shift of S 3i / S^f L = 1.04 and S^/S^ A AC = 0.95 using the 
purely Adelberger and NACRE rate compilations, respectively; the smallness of these shifts 
just restates the success of the SSM in light of the SNO observations. 

If we adopt the scaling laws and propagate the errors according to the usual rules, we 
have 

This gives a dispersion of c^/,!^ = 0.24 for both compilations. This is much larger than 
the small shifts in the mean found in the above paragraph. Moreover, we see that to solve 
the BBN 7 Li problem with reaction rate uncertainties alone requires a ~ 2a change in 5*34. 
Thus we find that this solution to the 7 Li problem is excluded at ~ 95% CL. We now turn 
to numerical results which will confirm and better quantify this limit. 
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4.2 Numerical Formalism and Results 



Our analytic discussion uses standard error propagation which is good only to first order, 
and assumes gaussian errors as well as linearity. To explore this scenario more rigorously, 
we perform this calculation numerically, taking into account the non-gaussian nuclear errors 
and non-linear scalings. We set out to perform a Monte Carlo integration of an integral of 
the form: 



where S is a set of reaction rates (such as the rates already listed) and <ps is the 8 B solar 
neutrino flux. Cssm, C-nuc, and Csno are the likelihood distributions of the Standard Solar 
Model given a reaction network and a solar neutrino flux, the reaction network given various 
rate compilations, and the total 8 B solar neutrino flux given by the SNO collaboration. 

In order to test the reliability and accuracy of this method, we first predict the total 8 B 
neutrino flux given a complete reaction network, using both the Adelberger and NACRE 
compilations and then compare to the predictions shown in the works of Bahcall et al [23J. 
The integral we are performing is: 



A Monte Carlo integration uses one of the likelihood functions to draw random numbers and 
average the remaining function over those generated random numbers. For our case, we will 
generate random numbers for the independent reaction rates given by either the Adelberger 
or NACRE compilations. We combine statistical and systematic uncertainties by adding 
them in quadrature, since only the total uncertainty is needed for this analysis. We generate 
gaussian or piecewise gaussian distributions for the reaction rates, depending on whether the 
quoted errors are symmetric or asymmetric about the most likely value. For each random 
draw of the reaction rates, we can calculate a solar neutrino flux, given the scalings shown 
in eqn. Once a large sample of 8 is created, we can calculate its likelihood distribution. 
To summarize: 

1. Cnuc(S) generates reaction rates randomly. 

2- £ssm(S, 8 ) enforces the scalings in eqn. 

3. The resulting sample of 08 is used to find £(0s). 

The normalization or best value and the errors are calculated separately. The flux val- 
ues for Adelberger and NACRE, as given in Tables 7 and 9 in j22] are the standard solar 




(11) 




(12) 
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model predictions for the neutrino fluxes, adopting each compilations best fit values, with- 
out marginalizing over the reaction network. The errors are then propagated separately, as 
described in j2S| using the scalings already mentioned. The scalings are valid in determining 
the uncertainties to 10%. Thus, we will adopt the scalings shown in eqn. normalized such 
that when a given compilation is used, we reproduce the values listed in 



-0.40 / c \ 0.81 / c \ 1-0 / c \ -1.0 

^ = 5 .05xl0«(^r) |^] Q|U ,13) 





-0.40 / c \ 0.81 / c \ 1-0 / c \ -i.o 
■ vw - 544 xl0 6 ( ° U \ ( ° 33 \ ( bu ) ( bl7 ) f-^L.) (U) 

<j.<±<± ^ iu ( oN AC I \ oNAC I oNAC \ oNAC I \ oADL I v 

\°34,0 / 

By using the Adelberger scaling relation to predict the NACRE scaling relation and vice 
versa, we can verify the accuracy of these fits. We find deviations from the relations listed 
above at the 8 or 9% level, thus we adopt an overall 10% systematic uncertainty in the pre- 
dicted flux. Also, since the resulting distributions are non-gaussian, we expect our marginal- 
ized best fit neutrino fluxes to be different from the neutrino flux determined by adopting 
only the best values of the reaction rates. 

We find remarkable agreement between our confidence intervals and those placed by 
Bahcall et al ■ Our results are summarized below, as well as in figure 



b$ DL = 5.09 
^ AC = 5.19 



1 -°-OM?oS){ stat ) ± O.lO(syst)} x 10 6 cm" 2 s" 1 (15) 



-0.16(0.29)^ 

l -°-lfi ( S( stat ) ± O.lO(s^)] x 10 6 cm" 2 s" 1 , (16) 



where the flux numbers listed are the most likely values for the Adelberger-based and 
NACRE-based compilations and their respective 68% (95%) confidence limits, as deter- 
mined from the marginalized likelihood distributions. Notice that our most likely values are 
different than the fluxes determined by adopting the best values for the reaction rates. This 
shift in best values is due to the marginalization over the non-linear scalings and asymmetric 
nuclear errors. Had the scalings been linear and additive, and all nuclear errors symmetric, 
no shift would have been seen. As one expected, NACRE has slightly inflated errors as 
compared to Adelberger. This is simply due to NACRE's overall larger rate uncertainties, 
as shown in our analytic work. 

As discussed earlier, the solar neutrino flux depends primarily on the S17 and S34 reac- 
tions. The S33 and reactions have little impact on the results due to their small errors 
and the weak dependence of the flux on them. The S11 rate has negligible effect in the 
Adelberger compilation, but has significant impact in the NACRE compilation's results. 
NACRE's uncertainty for this rate is a factor of 2 larger than the Adelberger's compilation. 
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Figure 1: Shown are the Standard Solar model predictions of the total 8 B neutrino flux. The 
binned likelihood based on the Adelberger (NACRE) rate compilation is plotted with solid 
(open) squares. The 10% systematic error has not been included here. 

Below we will use the differing results of these two compilations as a probe of the Su error 
assignment. 

Given the scalings in eqn. El we can use the SSM and the SNO measurement of the total 
8 B neutrino flux to constrain these rates in the following combination: x = S^S^f 1 . As 
before, we will generate random numbers for the independent reaction rates given by either 
the Adelberger or NACRE compilations. However, we now fix S17 and S 34 with various 
values of x. For each random draw of the reaction rates, we can calculate a solar neutrino 
flux, given the scalings shown in eqn. |S] With this flux, we then average Csno^s) over the 
sample to find the likelihood of a given x. To summarize: 

1. Cnuc(S) generates reaction rates randomly. 
2- CssAi(S,(ps) enforces the scalings in eqn. 
3. Calculate C(x) = (£sjvo(0s))- 

Using the Sn and S33 from the Adelberger and NACRE compilations respectively and the 
S e 7 from Adelberger compilation, and the SNO collaborations constraint on the total 8 B 
neutrino flux, we place the following constraints on x. 

[keV b] L81 (17) 



x 



ADL 



0.0119 



1.0 



+0.16 (0.35) 
-0.14 (0.25) 
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Figure 2: Shown are the likelihood distributions of the parameter x = SnS^f , given the 
subset of reactions from the Adelberger (solid) and NACRE (dashed) compilations respec- 
tively and the SNO collaborations measurement of the solar 8 B neutrino flux. The 10% 
systematic uncertainty in the flux scalings has been included as gaussian. 

[keV b] 1 - 81 , (18) 

where the most likely values, the 68% (95%) confidence intervals. The 10% systematic error 
has been included in the calculation and assumed to be gaussian. These resulting likelihoods 
for x are shown in figure 121 The Adelberger and NACRE -based results agree quite well 
with each other. With differences mainly attributable to the larger error in Su adopted by 
NACRE. 

Since we are constraining x only, we cannot determine the Sn and S34 reactions uniquely. 
We require additional information. If a total p — p or 7 Be neutrino flux measurement existed, 
we could in principle determine both cross sections. Since we are using the Sun to constrain 
systematic errors in the normalization of S34, in an attempt to fix the BBN 7 Li problem, we 
will adopt various experimentally-determined values of Sn to place constraints on S34. Once 
a value of Sn is adopted, we convolve the x likelihood distribution with the experimental 
Sn distribution to get our S34 likelihood. 

Besides using the Adelberger and NACRE rate compilations for 5*17, we also use two 
more recent determinations. We use the recommended values from Junghans et al. [39J, and 
Davids and Typel |%0] . The Junghans quoted value, Sn = 21.4 ± 0.5(expt) ± 0.6(theor) 
eV b, is based on several direct capture data sets. The Davids and Typel value, Sn = 



x 



NAC 



0.0121 



1.0 



+0.21 (0.46) 
-0.17(0.32) 
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Figure 3: Shown are the likelihood distributions of S34, given Sn measurements listed in 
Table EJ The upper (lower) panel shows the results using the Adelberger (NACRE) compi- 
lation for the Sn, S33 and S e 7 reactions. We have used values for Sn from Junghans [HH] 
(solid), Davids (dashed) and Adelberger |23J and NACRE j2H] (dotted). Again, the 10% 
systematic uncertainty in the scalings has been included and assumed gaussian. 
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Table 4: Shown are the constraints placed on S'34 using reaction rates from various sources. 
Column 1 lists the adopted Sn constraint used, while Columns 2 and 3 show the compilation 
used for the Sn and S 33 reaction rates. The S'34 numbers cited are the most likely values 
and their 68% (95%) confidence intervals. 



Adopted Sir (eV b) 


Adelberger-based 24 


NACRE-based 125 


Adelberger 24 
Sir = 19.012:8 


C _ n p- 1+0.15 (0.34) 
D34 — U.Ol_ i2( 21) 


N.A. 


NACRE [23 
Sir = 21.0 ±2.31 


N.A. 


C n n +0.17(0.38) 

034 — u.di_ 12(0 22) 


Junghans |29_] 

Sir = 21.4 ± 0.5(expt) ± 0.6(theor) 


C _ n /i o+O.lO (0.23) 
^34 - U - 4 »-0.08(0.15) 


q _ n ,0+0.14(0.30) 

^34 - u.4y_ 011(019) 


Davids |H| 

Sir = 18.6 ± QA(expt) ± l.l(extrp) 


C _ n p- 7 +0.13 (0.30) 
°34 — u -° '-0.11(0.19) 


C _ n ^0+0-17(0.39) 
°34 — u - oy -0.13(0.24) 



18.6±0A(expt) ± l.l(extrp) eV b, is based on both direct capture and Coulomb dissociation 
measurements, excluding the Junghans data set because it is systematically higher than the 
other data sets. Had the Junghans data been used, the value of S i7 would lie between the 
two cited values. We will adopt the cited numbers, keeping in mind that the difference in 
their values are a measure of this systematic difference. 

Our constraints in Tableware based on the likelihood functions in figure El We find that, 

S 34 > 0.35 keV barn (19) 

at 95% CL for the case of the NACRE Sir value. Other choices give slightly higher limits, 
e.g., Adelberger with the Davids Sir gives S 34 > 0.42 keV barn. 

As shown in Table El these limits on S34 place essentially identical limits to 7 Li production 
in BBN. Thus, eq. ()19|) . along with the fiducial BBN results in Tabled demands that 

C^f) > x ID" 10 , (20) 

\ n / BBN 

where we have fixed the reaction normalization such that S34 = 0.35 keV barn, but propa- 
gated the other nuclear uncertainties in the BBN code and convolved the predictions with 
the WMAP determination of the baryon density pQ . We see that this alleviates the BBN 7 Li 
problem somewhat, but still requires a combination of effects to fix the problem-i.e., that 7 Li 
observations be systematically low, in addition to adopting the limits to nuclear systematics 
we have derived. 
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Put another way, we can ask how far a "nuclear-only" BBN solution stretches our con- 
straints on 3 He(a,7) 7 Be systematics. We saw in §3.21 that for halo star observations to 
reflect the primordial 7 Li abundance requires that S34 be systematically lowered, to 53% 
and perhaps 27% of its fiducial value. A reduction of Sgf w < 0.267 keV barn is excluded 
at the 99.5% CL for the NACRE case (and above for others in Table HJ). A reduction of 
S3™ < 0.136 keV barn is excluded at more than 99.9999% CL. This restates our finding that 
the solar constraints on 5*34 remove this reaction as the main suspect in the 7 Li problem. 

5 Discussion and Conclusions 

The hot big bang cosmology has seen a great triumph in the agreement between the baryon 
density found by WMAP and the BBN value implied by the D/H ratio measured at high 
redshifts. However, this triumph is somewhat muted by the much poorer agreement between 
the primordial 7 Li value as predicted from BBN and the WMAP baryon density, and the 
observed values seen in halo stars. The predictions are at least a factor of 2 above the 
observations. This discrepancy impels a search for any possible systematic errors, which 
could either explain the mismatch, or if no systematics can be found, would reveal the true 
seriousness of the problem and a need for a more fundamental solution. 

In this paper we have considered the effect of systematic errors in the nuclear reactions. 
In particular, we have focused on the 3 He(a,7) 7 Be reaction, which is the sole important 
production channel of 7 Li at the WMAP baryon density. As such, systematic errors in this 
reaction have an immediate impact on the BBN 7 Li abundance. And indeed, while there has 
been extensive and careful work for this reaction, both fronts meet with technical difficulties 
which leave open the possibility for systematic errors in the absolute normalization of this 
rate. 

Thus we have identified a new constraint on this reaction, coming from its influence on 
7 Be and 8 B production in the Sun, and the associated 8 B solar neutrinos. The excellent 
agreement between the standard solar model and the total measured 8 B neutrino flux places 
demands that the underlying nuclear reactions cannot have large systematics. In particular, 
using the solar neutrino theory and observations, as well as some information on other 
reactions, notably 7 Be(p, 7) 8 B, we find that S34 cannot be smaller than 65% of its fiducial 
value (e.g., NACRE or Adelberger). This limit is strong enough to exclude the 3 He(a,7) 7 Be 
reaction as the dominant solution to the BBN 7 Li problem. 

Other nuclear solutions to the 7 Li problem are logically possible but in fact unlikely. 
While many reactions are important for 7 Li production, the requirements that we not spoil 
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agreement with D, and not (further) underproduce 4 He, leads us to focus on reactions which 
only affect 7 Li. Since we have shown that the production channel cannot be lowered suffi- 
ciently, we might hope to increase 7 Be destruction. This is done via the 7 Be(n,p) 7 Li reaction, 
followed by 7 Li(p, 7)2 4 He. The Sun does not constrain 7 Be(n, p) 7 Li because the solar interior 
has a negligible neutron density. However, this reaction is nevertheless very well-studied be- 
cause its inverse is a common laboratory neutron source. Since deuterium observations and 
CMB determinations suggest a baryon density on the high side, the destruction of 7 Li through 
the reaction 7 Li(p, 7)2 4 He has negligible impact. Its mirror reaction, 7 Be(n, 7)2 4 He, impor- 
tant on the higher baryon density side, is negligible compared to 7 Be(n,p) 7 Li. Furthermore, 
7 Li has a somewhat weaker dependence on the destruction cross section ( 7 L1bbn oc S^f 5 S~^ 74 
[01), so that the needed systematic error would be even larger than what we have considered 
for the production channel. 

Thus nuclear solutions do not seem allowed by the current data. Of course, it remains 
possible that extremely large (factors > 100) systematic errors lurk in otherwise negligi- 
ble 7 Li production and destruction channels |27j . For these reasons, continued efforts to 
improve nuclear cross section experiments and theory (with particular attention to abso- 
lute normalizations and systematics) will reap benefits for BBN as well as solar neutrinos. 
Tighter experimental errors (including systematics) will reduce the BBN theoretical uncer- 
tainty budget, which will not only further clarify the seriousness of the 7 Li problem, but 
also allow for stronger constraints on astrophysics ^T] when and if the 7 Li problem is re- 
solved. In this respect, we particularly call attention to the 3 He(a,7) 7 Be reaction, but also 
to 7 Be(jo, 7) 8 B, as they are undoubtedly linked through solar neutrinos. Determining a more 
accurate low-energy extrapolation in either of these reactions will impact the other through 
the solar neutrino constraint on the parameter x = S17S34 81 . 

Where, then, does the 7 Li problem stand? We have found nuclear reaction systematics are 
very unlikely to be the dominant source of the discrepancy. Of the remaining possibilities, the 
most conservative is that the problem is dominated by systematic errors in the observational 
7 Li value. This could either be due to difficulties in the understanding the stellar parameters 
and in extracting the abundance from spectral lines, or from stellar evolution effects which 
deplete Li without introducing large dispersion in the Spite plateau. A similarly conventional 
solution would ascribe the 7 Li discrepancy to a combination of nuclear and observational 
systematics, both at the edge of what is currently allowed. 

Finally, a more radical but intriguing possibility would be that new physics is required. 
If this is so, nature has been somewhat subtle in revealing this twist, as the perturbation to 
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standard BBN has been small enough not to be noticed until now. 3 Nonstandard scenarios 
have already been proposed to alleviate the 7 Li problem by introducing new physics, e.g., 
by a late-decaying gravitino [41 j . However, most of the scenarios require fine tuning, as one 
wishes to reduce 7 Li without spoiling the superb concordance between deuterium and the 
CMB. 

In summary, we use solar neutrinos to remove the possibility of a solution to the 7 Li 
problem from the 3 He(a,7) 7 Be reaction, and thereby cast more doubt that the problem 
is due to nuclear systematics. By removing a possible resolution, we have both clarified 
the problem, and made it more acute. In our view, the most important arena now is the 
observations and astrophysics which lead to the primordial 7 Li inference. And while we 
continue to suspect that this is the likely solution, a parallel examination of nonstandard 
BBN scenarios is at this point not unwise. 
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